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APPLICATION  OF  TWO-DIMENSIONAL  DISCRETE-ORDINATES 
METHOD  TO  MULTIPLE  SCATTERING  OF  LASER  RADIATION 


I.  INTRODUCTION 

Light  propagation  through  an  optically  thick  particulate  medium  is  basically  a 
multiple-scattering  problem  in  which  rays  or  photons  traverse  a  medium  of  scatterers  and 
undergo  many  scattering  events  before  escaping.  A  natural  framework  to  deal  with  this  type  of 
problem  is  provided  by  the  theory  of  radiative  transfer.  The  linear  Boltzmann  equation  -  in  the 
context  of  radiative  energy  also  termed  the  "equation  of  transfer"  -  governs  the  radiation  field 
in  a  medium  that  absorbs,  emits,  and  scatters  radiation. 

The  complexity  of  the  equation  of  transfer  forces  one  to  implement  numerical 
methods  of  solution.  Such  methods  seek  to  introduce  approximations  that  convert  the  integro- 
differential  form  of  the  equation  of  transfer  into  a  system  of  algebraic  equations  that  is 
amenable  to  solution  by  a  digital  computer.  The  most  direct  procedure  is  the  discrete- 
ordinates  approach  in  which  the  radiance  distribution  function  I(? ,\)  is  replaced  by  a  discrete 
set  of  values  at  a  discrete  set  of  phase  space  points  A^).  The  derivatives  and  integrals 

appearing  in  the  equation  of  transfer  must  also  be  replaced  by  a  corresponding  discrete 
representation  using  finite  difference  and  numerical  integration  schemes.  In  this  way  one 
arrives  at  a  set  of  algebraic  equations  that  can  be  solved  by  using  standard  numerical  algorithms 
on  a  digital  computer. 

Unfortunately,  such  a  calculation  becomes  an  immense  undertaking  if  only  a  "brute 
force"  discretization  of  the  equation  of  transfer  is  employed.  For  example,  a  typical  mesh  size 
of  100  by  100  by  100  space  points,  10  wavelength  groups,  and  10  discrete  angles  would  yield  an 
extremely  large  set  of  simultaneous  algebraic  equations  involving  at  least  10°  discrete  phase 
space  points,  which  is  a  rather  formidable  task  even  on  the  CDC-7600  or  Cray-1  computers. 

The  discrete-ordinates  method  is  discussed  in  detail  in  references  1-5.  In  the 
context  of  neutron  transport  theory,  a  comprehensive  treatment  is  contained  in  the  monograph 
of  Bell  and  Glasstone.  Modem  neutron  transport  methods  have  only  recently  been  applied  to 
atmospheric  radiative  transfer  problems  by  Gerstl. '  Following  reference  5,  we  summarize  the 
main  features  of  the  method: 

(a)  The  solution  of  the  transfer  equation  may  be  derived  explicitly  and  therefore 
the  intensity  and  flux  computations  do  not  only  depend  on  the  total  optical  depth  of  aerosol 
layers. 


(b)  The  method  yields  the  radiation  field  at  all  spatial  locations  as  well  as  the 
reflection  and  transmission  without  additional  computational  effort. 

(c)  The  computer  time  required  to  calculate  intensity  and  flux  is  relatively  small  as 
compared  with  other  techniques. 

For  two-  and  three-dimensional  geometries,  the  solution  of  the  discrete-ordinates  equations 
becomes  increasingly  complex.  Extensions  from  two  to  three  dimensions  result  in  severe 
computational  problems  because  of  the  sheer  size  of  the  system  of  equations  rather  than  any 
basic  complications  in  numerical  solution  algorithms  from  the  added  dimension.  A  discrete- 
ordinates  code  THREETRAN8  which  is  capable  of  treating  three-dimensional  geometry,  exists 
only  in  an  experimental  version. 


When  looking  for  rigorous  solutions  to  the  problem  of  light  scattering  amt 
propagation  in  dense  aerosols,  it  thus  becomes  imperative  to  explore  the  usefulness  and  the 
range  of  applicability  of  the  existing  two-dimensional  transport  codes.  We  note  that  radiative 
transfer  in  two  dimensions  has  earlier  been  studied  by  Fowler^  and  Fpwler  and  Sung.  ®  The 
purpose  of  this  paper  is  to  apply  the  two-dimensional  code  TWOTRAN**  to  the  description  of  a 
laser  beam  scattered  (reflected)  from  a  fixed  target.  This  is  obviously  relevant  to  the 
evaluation  of  the  effectiveness  of  electro-optical  systems  under  aerosol-loaded  and  adverse 
weather  conditions.  After  summarizing  the  main  features  of  the  TWOTRAN  code  and  the 
solution  algorithm  in  section  2,  we  make  use  in  section  3  of  the  advantages  of  the  delta  M 
method.  In  section  4,  we  describe  our  model  scenario  where  a  collimated  laser  beam 
propagating  in  a  homogeneous  rural  aerosol  is  subsequently  scattered  from  an  isotropic  target. 
The  average  intensity  and  Beer/Lambert  correction  factor  are  introduced  in  section  5.  They 
allow  us  to  make  a  comparative  study  of  results  in  the  rectangular  (x-y)  and  cylindrical  (r-z) 
geometries,  which  is  contained  in  section  6.  Section  7  is  devoted  to  the  discussion  of  the  results 
obtained,  and  conclusions  outlining  future  lines  of  research. 


2.  SOLUTION  ALGORITHM 

For  calculations  reported  here,  the  general-purpose  code  TWOTRAN-II  has  been 
adapted  to  radiative  transfer  calculations  of  a  laser  beam  in  a  realistic  aerosol  medium. 

TWOTRAN-n  is  a  general  purpose  FORTRAN  code  that  solves  the  two-dimensional 
multigroup  transport  equation  in  (x,y),  (r,Ol,  and  (r,z)  geometries.  Both  regular  and  adjoint, 
inhomogeneous  and  homogeneous  problems  subject  to  vacuum,  reflective,  periodic,  white  or 
input-specified  boundary  flux  conditions  are  solved.  General  anisotropic  scattering  is  allowed 
and  anisotropic  inhomogeneous  sources  are  permitted. 

The  discrete  ordinates  approximation  for  the  angular  variable  is  used  in  finite 
difference  form  which  is  solved  with  the  central  (diamond)  difference  approximation.  Negative 
fluxes  are  eliminated  by  a  local  set-to-zero  and  correct  algorithm.  Standard  inner  and  outer 
iterative  cycles  are  accelerated  by  coarse  mesh  which  may  be  independent  of  the  material 
mesh,  .pur  computations  with  TWOTRAN  performed  with  a  pointwise  convergence  criterion  of 
2  x  10  6  using  144  discrete  directions  (S j^)  and  42  by  42  mesh  point  required  the  running  time  of 
1.2-1. 5  min  on  the  CDC-7600  machine. 

In  the  absence  of  volumetric  sources,  the  equation  of  transfer  to  be  solved  is 

Q-Vl(r.n)  +  ot(r)I(r,0)  =  /  os(r,p)I<r,0' )dD'  {1) 

where 

a  t  is  the  volume  extinction  coefficient  and  I(r,  S2 )  is  the  radiance  distribution 
function.  In  nuclear  radiation  transport  terminology,  <xs(r,  p)  is  the  differential-scattering- 
transfer  cross  section,  where p  is  the  cosine  of  the  angle  between  the  directions  Q' and  fi  of  the 
incoming  and  outgoing  pencils  of  radiation.  In  terms  of  vg(r,p),  the  scattering  phase  function, 
P(r,p),  is  given  by 

P(r.p)  =  os(r,p)/oJ(r)  (2) 

with  the  normalizaton 


J  P(I.P)$  =  1  (3) 

J  s(r,p)dp  =  o°(r)  (4) 
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where 


i 


i 


,7  >s(r)  is  the  volume  scattering  coefficient  at  location  r. 

In  t  w  o-d  i  m  e  nsi  o  an  I  x-y  ami  i—  z  geometries,  the  TWOTRAN  code  requires  the  sets  of 
spatial  meshes  1  or  (r._,z  n)  and  discrete  directions  {  Qm|.  We  always  choose  the  ben  n  axis 

in  the  y  or  z  directions,  in  x-y  and  r-z  geometries,  respectively.  As  the  set  o'  b- 
directions  |  is  determined  by  the  choice  of  a  Gaussian  quadrature,  there  will  be  usually  no 
direction  available  along  the  y  or  z  axes.  It  is  therefore  convenient  to  divide  the  total  radi  <  i  :■ 
into  two  parts,  the  reduced  I  •  and  the  diffuse  radiance  Tj.  The  reduced  radiance  satisfies 
equation  (1)  with  the  right-hand  side  set  equal  to  zero.  In  (he  x-y  geometry,  for  example,  it  is 

lr.(0,r)  =  1(0, x  *  y0x/0y,0)  exp[-  ^  J  °t(*  -  y'0x/0y,y' )dy' ] 


w  ior 


TV-’-'  -  *  >  "y  o 

ft  and  SI  are  the  components  of  the  vector  ft  in  the  plane  x-y. 
x  y 

The  diffuse  radiance  satisfies  the  equation 

n*VId(r,0)  +  ot(r)Id(r,0)  =  J  os(r,p)Id(r,0' )dft'  +  Q(r,Q) 


where  Q  is  the  sour  :  *  function 


«(r.Q)  *  J  os(r.M)iI.i(r.n')‘1!1' 


(51 


(51 


(7) 


I 


I 


generated  by  the  reduced  radiance. 

In  order  to  obtain  the  solution  for  the  diffuse  radiance  from  the  TWOTRAN  program, 
both<rs  and  Q  should  be  represented  by  finite  Legendre  polynomial  expansions.  The  function 

*s(r,M)  iS  vritten  as 


osd,P)  =  .Io  oJCDP^P) 
d 


(«) 


where  the  expansi on  coefficients  are  assumed  to  be  known.  To  write  do. 7  i  .u  \xpansion  for 
the  source  function  O,  we  take  the  case  of  rectangular  x-y  geometry  as  an  example.  For  a 
colli  nated  laser  beam  with  the  spatial  shape  described  by  a  function  F(x)  at  y  =  0 


I(x,y=0,Q)  =  F(x)6(n  -  y) 


(9) 


by  virtue  of  equation  (5)  the  reduced  radiance  becomes 


I  .(ft,r)  =  F(x)S(ft  -  y)exp[-  J  ot(x,y‘)dy'] 
n  -  -  o 


Equation  (7)  now  yields 


Q(r,n)  =  F(x)c  (r,n-y)exp[-  /  at(x,y')dy'] 


do) 


(in 


V*  is  evident  bv  inspection  of  equations  (8)  and  (11),  the  expansion  of  tlie  source 
function  involves  the  series  expansion  of  equation  (11)  with  th  ;  Legendre  polynomials  Pj(ftft) 
taken  at  the  argument  P n  ( ).  For  strongly  asymmetric  phase  functions,  the  computational 
dif ficulties  arise  because  such  phase  functions  cannot  be  represented  by  polynomials  of  low 
degree.  A  way  out  of  the  predicament  is  the  delta  M  method,  which  we  will  now  discuss  in  tee 
context  of  the  TWOTRAN  transport  code. 
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3. 


THE  DELTA.  M  METHOD  FOR  THE  TWOTRAN  CODE 


The  delta  M  method,  which  relies  essentially  on  matching  the  first  2M  phase 
function  moments  and  using  a  Dirac  delta  function  representation  of  forward  scattering,  is  an 
extension  of  the  6  -  Eddington  approximation.  '  The  method  has  recently  been  discussed  i  i 
connection  with  the  scaling  group  of  the  radiative  transfer  equation.1  In  order  to  adapt  the 
delta  M  method  for  our  purposes,  we  note  that  the  only  physical  quantities  required  to  run  the 
TWOTRAN  code  are  the  expansion  coefficients  0  ^  s  in  equation  (8)  and  the  extinction 
coefficient  at  each  mesh  side.  When  the  phase  function  is  approximated  by  a  Dirac  delta 
function  forward  peak  and  a  (2M  -  l)-term  expansion,  the  function  0^  in  equation  (8)  becomes 


*os(r,M) 


2M-1 

I 

£=0 


21 


+  1*0 

ir  vr)P,(M) 


(12) 


where  for  highly  asymmetric  phase  functions  the  number  2M-1  will  be  much  smaller  than  L. 
The  condition  that  the  first  (2M-l)-phase-function  moments  be  matched  exactly  leads  to 


with 


f 


°sM(r> 


Since  the  extinction  coefficient  transforms  as 


f 


(13) 


(14) 


(15) 


the  new  equation  of  transfer  for  the  transformed  radiance  I*  becomes 


Q-VI*  +  *0tI*  =  J  *as(r,p)I*(r,n,)dQ'  (lb) 

Thus,  by  comparing  equation  (16)  with  equation  (1),  we  see  that  the  delta  M  method  does  not 
change  the  functional  form  of  the  transfer  equation. 

Actually,  as  explained  in  section  2,  we  solve  the  equation  of  transfer  for  the  diffuse 
radiance.  This  is  obtained  from  equation  (6),  in  which  the  og  and  0  t  have  been  transformed 
according  to  equations  (12)  and  (15).  The  reduced  part  of  the  radiance  is  still  given  by  a  simple 
analytic  expression,  obtained  by  combining  equations  (5)  and  (15): 

In*  =  fT  *  <foxKx  *  )dy' 3  (17) 

17 or  a  homogeneous  medium,  with  the  collimated  laser  ’to  a  propagating  along  the  y 
axis,  equation  (17)  simpifies  to 

Iri  =  Iri  exp(fo°y)  (18) 

Therefore,  from  the  approximate  equality 


*  * 

Ir4  +  S  Iri  + 
n  a  n  d 


(19) 
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c  =  J-  =  exp(fojy)  -  1  +  yfi-  (20) 

Ari  Vi 

Equation  (20)  enables  one  to  evaluate  the  multiple-scattering  corrections  to  the  exponential 
decay  expressed  by  the  Beer-Lambert  law. 

4.  MODEL  SCENARIO 

In  this  section  we  describe  the  results  pertaining  to  a  simple  scenario  in  which  the 
scattering  medium  is  formed  by  a  homogeneous  distribution  of  99%  relative  hujnidity  rural 
aerosol.  The  rural  model,  as  developed  by  the  Air  Force  Geophysics  Laboratory,  5  is  intended 
to  represent  the  aerosol  conditions  one  finds  in  continental  areas  which  are  not  directly 
influenced  by  industrial  aerosol  sources.  The  optical  characteristics  of  the  medium  are 
generated  with  the  aid  of  the  AGAUSX  code,  as  described  in  references  16  and  17.  Figure  A-l 
(appendix  A)  shows  relative  locations  of  the  laser  source  and  the  isotropically  reflecting  target 
in  the  x-y  (r-z)  plane.  The  x  (r)  axis,  named  also  "height",  is  chosen  perpendicular  to  the  surface 
of  the  earth,  while  the  y  (z)  axis,  named  "distance"  is  taken  along  the  axis  of  symmetry  of  the 
laser  beam.  In  r-z  geometry,  the  medium  has  an  azimuthal  (rotational)  symmetry  around  the  z 
axis.  The  symmetry  condition  in  x-y  geometry  may  be  illustrated  by  imagining  that  the  medium 
is  infinite  in  extent  in  the  direction  perpendicular  to  the  x-y  plane.  Both  the  x  (r)  and  y  (z) 
directions  are  divided  into  six  coarse  mesh  cells  whose  boundaries  are  given  by  x  (r)  =  0.0,  0.002, 
0.010,  0.050,  0.8,  1.0,  2.5  kra,  and  y  (z)  =  0.0,  3.0,  3.8,  4.0,  4.006,  4.106  and  4.4  km.  In  addition, 
each  coarse  mesh  cell  is  subdivided  into  fine  mesh  cells  of  equal  width.  In  the  x  (r)  direction, 
the  sequence  2,  8,  10,  10,  2,  10  specifies  the  number  of  fine  meshes  per  coarse  mesh.  The 
corresponding  numbers  in  the  y  (z)  directions  are  20,  8,  5,  2,  4  and  3.  Thus  the  x-y  (r-z)  plane 
has  been  covered  by  42  x  42  fine  mesh  cells  with  unequal  sides.  In  x-y  geometry,  the 
coordinates  of  the  laser  source  are  xQ  =  1.5  m  and  yQ  =  3.0  km.  The  source  emits  a  Gaussian 
collimated  beam  with  the  radiance  distribution  for  y  =  yQ  given  by  the  formula 


I(x,y0,fi)  =  F/(n1/2w)exp[-(x-x0)2/w2]6(ft-y)  t 


where  F  is  the  total  energy  of  the  beam,  and  5  denotes  the  Dirac  delta  function.  In  realistic 
situations,  one  usually  deals  with  a  pulse  source.  Typically,  the  emitted  pulses  are  of  1  M*.v 
power  and  20  nsec  duration.  We  thus  set  F  =  2  x  10~ZJ  per  pulse  and  take  w  =  0.5  m.  As  the 
distance  between  the  source  and  the  target  is  equal  to  1  km,  a  point  source  having  the 
divergence  of  1  mrad  would  produce  a  spot  of  1  m  in  diameter  at  the  target's  location.  The 
collimated  Gaussian  beam  with  half  width  w  =  0.5  can  therefore  be  viewed  as  a  beam  equivalent 
to  that  produced  by  a  divergent  point  source,  at  least  as  far  as  the  target  illumination  is 
concerned.  The  target,  with  assumed  reflectivity  r.  =  0.5,  reradiates  the  incident  beam 
isotropically  generating  a  diffuse  radiation  field  in  tne  x-y  plane.  The  ground  reflectivity 
(surface  albedo)  for  x  -  0  is  assumed  to  be  r  =  0.1.  In  r-z  geometry,  the  radiance  distribution 
function  in  the  plane  z  =  zQ  is  given  as  b 

I(r,z0,ft)  s  exp(-r2/w2)6(ft  -  z)  (22) 

nw 

Equation  (22)  implies  that  the  beam  axis  coincides  with  the  z-axis.  As  before  we  take  zQ  =  3  km 
and  suppress  the  ground  reflectivity  r  .  For  rectangular  x-y  geometry,  as  the  beam  is  parallel 
to  the  ground  surface  x  =  0,  the  exact® value  of  r  has  only  a  negligible  effect  on  the  results  of 
our  calculations. 


..  . 


We  choose  the  concentration  of  the  aerosol  particles  to  be  described  by  the  number 
density  of  50,000  particles  per  cm  .  This  corresponds  to  the  extinction  coefficient  O  =  1.32 
km*  at  the  wavelength  A  =  1.06  JJIB  .  For  A  -  0.55  pH*  the  extinction  coefficient  is  2.57  km 
The  meterological  (visual)  range  V,  as  defined  by  the  Koschmieder  formula 


V  ~  o(A=0.55  m^)  ’  (23) 

'l 

is  then  V  =  1.52  km.  As  the  number  density  increases  to  100,000  particles  per  cm'  ,  the 
extinction  coefficient  becomes  0=  2.64  km*  and  the  meteorological  range  decreases  to  V  =  0.76 
km.  The  optical  depths  at  A  =  1.06  pm  for  the  distance  between  laser  and  the  target  (1  km)  in 
the  two  cases  considered  are  1.32  (for  V  =  1.52  km)  and  2.64  (for  V  =  0.76  km). 


5.  AVERAGE  INTENSrrY  AND  BEER-LAMBERT  CORRECTION  FACTOR 

The  knowlege  of  the  radiance  distribution  function  l(r, fallows  us  to  introduce  the 
average  intensity 


<I(r)>  =  ^  J  I(r,3)d0 


(24) 


The  average  intensity  does  not,  in  general,  represent  the  power  flow  but  is  proportional  to  the 
radiant  energy  at  a  spatial  mesh  point  r.  In  equation  (24)  and  in  the  remainder  of  this  work, 
I(r,i$)is  the  diffuse  radiance  at  a  two-dimensional  spatial  mesh  point  specified  by  the  radius 
vector  r.  The  direct  (reduced)  beam  intensity  is  omitted. 


In  order  to  investigate  the  relative  importance  of  multiple-scattering  corrections  to 
the  Beer-Lambert  law,  we  introduce  a  correction  factor  C  which  is  entirely  due  to  multiple 
scattering.  As  the  target  reflects  isotropically,  the  reduced  (unscattered)  light  intensity  at  a 
distance  from  the  target  is 


I 


r 


V*1 

*Z7 


(25) 


where  PQ  is  the  power  emitted  by  the  target  and  I  =  o£  the  optical  depth  between  the  target 
and  a  receiver.  The  total  intensity,  a  sum  of  the  reduced  and  diffuse  parts, 


I  =  Ir  +  Id 


(26) 


may  be  written  in  terms  of  a  correction  factor  C  to  Beer-Lambert's  law  as 


P_6 
I  « 


-I 


4  nl 


ip  (1  +  C)  ,  so  that 


Solving  equation  (27)  for  C,  we  obtain  explicitly 


(27) 


I 


c  = 


d  _t 

— r  e 


P0/(47i£‘) 


(28) 


As  the  diffuse  intensity  is  created  within  the  medium  due  to  multiple  scattering,  the  correction 
factor  vanishes  when  only  target  reflection  and  absorption  are  considered. 
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6. 


NUMERICAL  RESULTS 


The  quantities  of  interest  introduced  in  section  4,  namely  the  average  intensity, 
scattered  radiance  and  the  B/L  correction  factor  will  now  be  computed  in  x-y  and  r-z 
geometries. 

6.1  Rectangular  Geometry. 

Figure  A-2  shows  the  average  intensity  distribution  for  multiply-scattered  laser 
radiation  corresponding  to  the  visual  range  V  =  1.52  km.  In  figure  A-3,  the  same  result  is  shown 
with  the  spatial  scale  distorted.  This  displays  detailed  intensity  variations.  The  distortion 
represents  every  (nonequidistant)  spatial  fine  mesh,  as  defined  in  figure  1  and  section  4,  as  if  it 
were  equally  spaced.  This  feature  causes  the  intensity  peaks  created  by  the  aerosol  environ¬ 
ment  and  the  target  to  be  broadened  and  more  easily  recognizable  in  figure  A-3. 

In  figures  A-4  and  A-5,  the  average  intensity  distribution  is  shown  for  a  situation 
corresponding  to  V  =  0.76  km  where  the  aerosol  concentration  is  increased  to  100,000  particles 
per  cm  .  The  bimodal  distribution  in  figures  A-2  and  A-3,  with  the  target  radiation  being  the 
dominant  feature,  becomes  now  a  monomodal  and  the  target  is  no  longer  easily  recognizable. 
This  result  shows  that  for  sufficiently  high  aerosol  concentration  the  dominant  scattered  signal 
corresponds  essentially  to  the  interaction  of  laser  radiation  with  the  aerosol  environment. 

Dropping  thus  the  case  of  higher  aerosol  concentration,  we  list  in  tables  B-l  and  B-2 
(appendix  B)  the  B/L  correction  factor  for  V  =  1.52  km  assuming  a  detector  moving  in  a  plane 
parallel  to  the  surface  of  the  earth  and  pointing  toward  the  target.  In  actual  computations  the 
reduced  intensity,  as  defined  by  equation  (5),  was  obtained  by  executing  our  computer  code  in 
the  case  involving  no  scattering,  i.e.,  setting  the  scattering  coefficient  of  aerosols  v  =  0.  The 
coordinates  of  the  detector  were  determined  by  finding  a  ray,  in  the  set  of  discrete  directions, 
which  emanates  from  the  target  to  cross  the  line  along  which  the  detector  is  assumed  to  move. 

The  use  of  the  delta  M  method  reduces  significantly  the  computer  time  required  to 
obtain  the  radiance  distribution  function  and  the  correction  factors.  When  executing  the 
TWOTRAN  code  with  the  P^  expansion  (L=  5  in  equation  (8)],  the  CPU  time  was  90  sec  for  1764 
spatial  meshes  and  144  directions.  Employing  the  delta  4  approximation  results  in  the  reduction 
of  CPU  time  down  to  45  sec.  The  two  columns  of  the  correction  factors  in  tables  3-1  and  B-2 
refer  to  these  two  ways  of  computation. 

6.2  Cylindrical  Geometry. 

Leaving  unchanged  other  parameters,  we  proceed  to  the  cylindrical  (r-z)  geometry 
model.  Figures  A-6  and  A-7  show  the  average  intensity  when  the  visual  range  is  1.52  km,  while 
figures  A-8  and  A-9  correspond  to  the  visual  range  of  0.76  km.  As  before,  the  target  signal 
becomes  undistinguishable  from  the  aerosol-scattered  radiation  as  the  visual  range  becomes 
smaller  than  the  distance  between  laser  source  and  target.  Tables  B-3  and  B-4  list  the  scat¬ 
tered  radiance  and  the  B/L  correction  factors  for  the  aerosol  model  with  visual  range  V  =  1.52 
km  in  the  case  of  a  detector  moving  along  a  line  parallel  to  the  beam  axis. 


7. 


DISCUSSION 


Comparing  the  results  contained  in  the  figures  A-2  to  A-5  (x-y  model)  with  those  of 
figures  A-6  to  A-9  (r-z  model),  we  find  the  same  characteristic  features  in  the  diffuse  average 
intensity.  In  particular,  the  scattered  radiance  from  the  target  becomes  negligibly  small  as  the 
visual  range  is  shorter  than  the  distance  between  the  source  and  the  target.  'Ve  thus  conclude 
that  the  x-y  and  r-z  geometry  models  are  qualitatively  equivalent  as  far  as  the  target  detection 
is  concerned.  The  infinite  z  direction  in  x-y  geometry  is  only  an  expression  of  a  symmetry 
condition  and  is  not  a  perturbing  factor  as  long  as  results  are  considered  only  in  the  x-y  plane. 

To  the  extent  that  the  plots  of  average  intensity  show  a  qualitative  agreement 
between  the  two  geometries,  insight  into  the  quantitative  results  can  be  gained  by  inspecting 
tables  B-l  through  B-4.  For  smaller  optical  depths,  up  to  about  T  =  1.0,  the  agreement  in  the 
calculated  correction  factors  is  indeed  quite  satisfactory  and  suggests  a  functional  form  of  the 
correction  factor  described  by  a  universal  function  of  only.  As  the  distance  between  the 
target  and  the  detector  increases,  the  observed  discrepancies  can  be  attributed  to  the  ray 
effect.  This  phenomenon,  due  to  the  basic  discrete  ordinate  approximation  itself,  consists,  in 
essence,  in  solving  the  transport  equation  along  a  limited  number  of  discrete  characteristics 
(i.e.,  rays).  For  the  sake  of  illustration,  let  us  consider  an  isotropic  line  source  in  a  purely 
absorbing  medium.  Then,  clearly,  the  resultant  analytical  radiance  distribution  function  will 
possess  azimuthal  symmetry  about  the  source  line.  But  if  we  apply  the  discrete  ordinates 
method  to  this  situation,  the  resultant  radiance  will  consist  of  a  stepfunction  in  the  azimuthal 
angle  because  only  specific  directions  emanating  from  the  line  source  will  contain  the  source 
photons.  Thus  the  discrete  ordinates  method  approximates  the  azimuthally  uniform  radiance  by 
a  discrete  set  of  values  at  discrete  azimuthal  angles.  The  existing  remedies  like  those  trans¬ 
forming  the  discrete  ordinates  equations  to  spherical-harmonic  form,*®-^  will  be  an  object  of 
future  research. 


In  summary,  we  have  applied  the  discrete  ordinates  code  TWOTRAN  to  the  problem 
of  laser  beam  propagation  and  multiple  scattering  in  dense  aerosols.  The  delta  M  approximation 
has  proven  to  be  of  definite  advantage  for  our  applications.  The  main  result  of  our  investigation 
consists  in  validating  the  use  of  an  x-y  geometry  model  in  describing  a  realistic  three-dimen¬ 
sional  scenario.  This  is  achieved  by  computing  the  average  intensity  and  B/L  correction  factor 
in  x-y  and  r-z  geometries.  The  effect  of  aerosol  concentration  on  the  detectability  of  a  target 
follows  independently  of  the  choice  of  x-y  or  r-z  geometry.  In  addition,  as  we  have  shown  by 
evaluating  the  multiple-scattering  corrections  to  the  Beer-Lambert  law,  the  multiple-scattering 
contribution  to  the  target  signature  can  be  quite  dramatic  even  for  a  moderate  optical  thickness 
of  the  order  of  2.  As  the  optical  thickness  increases  beyond  T  =  1.5,  the  correction  factor  grows 
very  rapidly.  The  problem,  however,  in  the  context  of  the  discussed  scenario  becomes  less 
realistic  due  to  the  decreased  visual  range. 


Of  definite  interest  would  be  a  further  study  of  the  ray  effect  as 
comparison  of  the  results  obtained  with  the  results  employing  the  small  angle 
diffusion  approximations.  4 
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Rural  Aerosol 
Wavelength,  A=  !.08/xm 
Visual  Range  =  1.52km 
Max.  Intensity  =  3.32S-09 


Wavelength,  X-  1.06//m 
Visual  Range  =  0.76km 
Max.  Intensity  =  5.29E-10 


til  w. 


Figure  A-5.  Normalized  Average  Intensity  on  a  Distorted  Lattice  of 
Spatial  Mesh  Points  in  x-y  Geometry.  Visual  Range  0.76  km. 
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Rural  Aerosol 
Wavelength,.  A=  1.06/zm 
Visual  Range  =  1.52km 
Max.  Intensity  =  1.05E-06 


Rural  Aerosol 
Wavelength,  A=  1.06/zm 
Visual  Range  =  1.52km 
Max.  Intensity  =  1.05E-06 


Rural  Aerosol 
Wave  length,  A=  1.05/zm 
Visual  Range  =  0.76km 
Max.  Intensity  =  3.75Z-07 


Rural  Aerosol 
Wavelength,  X=  1.06/zm 
Visual  Range  =  0.76km  _ 
Max.  Intensity  =  3.78E-07 


Figure  A-9.  Normalized  Average  Intensity  on  a  Distorted  Lattice  of 
Spatial  Mesh  Points  in  r-z  Geometry.  Visual  Range  0.76  km. 
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APPENDIX  B 


TABLES 


Table  B-l.  X-Y  Geometry:  Target  Coordinates  x  =  0.001  km, 
y  =  4.000  km:  Detector  Height  =  0.200  km 


Detector's 
y  coordinate 


2.515 

3.529 

3.767 

3.822 

3.867 


Optical 

depth 


Radiance 


2.9792  x  10'1 
1.3051  x  10“} 
1.8615  x  10“} 
3.0306  x  10~: 
3.8526  x  10_1 


Beer/Lambert 
correction  factor 
P5  I  Delta  -4 


Table  B-2.  X-Y  Geometry:  Target  Coordinates  x  =  0.001  km, 
y  =  4.000  km:  Detector  Height  =  0.500  km 


Detector's 
y  coordinate 


Optical 

depth 


Radiance 


Beer/Lambert 
correction  factor 


P5  Delta  -4 


2.824 

3.556 

3.668 

3.904 


2.9857  x  10"2 
1.1178  x  10'1 
1.0735  x  10"J 
1.9973  x  10_l 


Table  B-3.  R-Z  Geometry:  Target  Coordinates  r  =  0.0  km, 
z  =  4.000  km:  Detector  Height  r  =  0.200  km 


Detector's 
y  coordinate 


Optical 

depth 


Radiance 


Beer/Lambert 
correction  factor 


P5 


2.959 

3.699 

3.775 

3.829 


1.398  8.9695  x  10“2 

0.477  5.1730  x  10"1 

0.397  9.4538  x  10'1 

0.348  2.8858  x  10+0 
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Detector's 
y  coordinate 


Optical 

depth 


Radiance 


Beer/Lambert 
correction  factor 
~P5  I  Delta  - 


